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ABSTRACT 

Based on 58 SLACS strong-lens early-type galaxies with direct total-mass and stellar-velocity dispersion 
measurements, we find that inside one effective radius massive elliptical galaxies with Mgff > 3 ■ 10'" Mq 
are well-approximated by a power-law ellipsoid with an average logaritmic density slope of (Tl^) = 
-d\og{pxot) I d\Qg{r) - 2.085;^QQj^ (random error on mean) for isotropic orbits with y6,- - 0, +0.1 (syst.) and 
CTy' < 0.20^Q Q2 intrinsic scatter (all errors indicate the 68% CL). We find no correlation of y'^^ with galaxy 
mass (Meif), rescaled radius (i.e. ^einst/^eif) or redshift, despite intrinsic differences in density-slope between 
galaxies. Based on scaling relations, the average logarithmic density slope can be derived in an alternative 
manner, fully independent from dynamics, yielding {y'^^ - 1 -959 + 0.077. Agreement between the two values 
is reached for (fir) - 0.45 + 0.25, consistent with mild radial anisotropy. This agreement supports the robust- 
ness of our results, despite the increase in mass-to-light ratio with total galaxy mass: Meu oc L^^^^^-^^^ . We 
conclude that massive early-type galaxies are structurally close-to homologous with close-to isothermal total 
density profiles (<10% intrinsic scatter) and have at most some mild radial anisotropy. Our results provide new 
observational limits on galaxy formation and evolution scenarios, covering four Gyr look-back time. 
Subject headings: gravitational lensing — galaxies: structure 



1. INTRODUCTION 

Understanding the internal structure of massive early- 
type galaxies (ETG) is essential, if we ever h ope to fully 
understand hierarchical galaxy formation (e .g. iDavis et all 
[1981 iFrenk et al.llT98l IWhite & FrenkHTggT ) and the com- 
plex interplay between dark-matter and baryons. A num- 
ber of tig ht relations, such as the Fundamental Plane (FP 



hereafter; iDressler et all Il987t iDiorgo vski & Davis' T987), 
the c o lor-m agnitud e relation ( e .g. [\^v anathan & Sandage 
119771; [Sandage & VisvanathanI 119781; iBoweretal.l 119921) 
and the relation betwee n black-hole and stell ar spheroid 
masses ([Kor mendv & Richsto nell995l;lM agorrian et aLlll 99*81 
lFerrarese & Merritti2000: .Gebhardt et al.i.200 0). indicate that 
there must be physical processes that are not dominated by 
stochastic processes of hierarchical galaxy formation. 

In the FP relation of ETGs both stars and dark-matter con- 
tribute to the structure and dynamics (cr^) and total mass-to- 
light ratio {MIL) . Hence, without knowledge of the dark- 
matter distribution and stellar populations, we remain clueless 
about whether changes in the FP are due to non-homology 
or due to changes in their stel lar mass- to-light, or both (e.g. 



Hiorth & Madsen 1995; Prugniel & Simien 1997;iTreu et al 



1999,: Gerhard et al. 2001; Berlin et al."2002; Cappellari et al 
200a [Nipoti et al. 2006; Graves et al. 2009). Given only first 



and second moments of the stellar velocity distribution, dis- 
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entangling these effects has remained difficult, because of 
the mass-anisotropy degeneracy: steepening of the density 
profile and changes in the orbital anisotropy can offset each 
other to yield similar kinematic profiles (e. g. lGerhard|[l993l: 
iGerhard et al.l ll998VLokas & Mamon"2003'). In studying the 
tilt of the FP (e.g. Renzini & Ciotti 1993), one can therefore 
not easily disentangle effects of mass structure (e.g. non- 
homology) from changes in stellar mass-to-light ratios, using 
only stellar kinematics^ 

^In recent papers (Bolton et al.' '2006"; 'Treu et al.' '20061 
Koopmans et al. 2006; Bolton et al. 2007; Gavazzi et al] 
20071: iBolton et al. 2008a), we analyzed a subsample of well- 



selected gravitational lens systems from the SLACS Sur- 
vey^, showing that massive elliptical galaxies have on av- 
erage close-to isothermal density profiles, with some mi- 
nor, but noticeable, intrinsic scatter between their logarithmic 
densi ty slopes (e.g. Koopmans et al. 2006; Treu & Koopmans] 
12004 . They follow the classical FP (Treu etal 2006), as 
well a s a tight Mass Fundamental Plane (MFP; Is'olton et all 
120071 r2008 b). where galaxy surface brightness is replaced 
by surface density. In all observable re spects, they foUovy 
the trends of normal elliptica l galaxies (Bolton et al. 2006t 
iTreu et all [2006; B olton et all 2008a; Tre u et al. 2QQ9) and 
these lens-based results can thus be exte nded to non-lens 
galaxi es in the same para meter space (e.g. iHvde & Bemardil 
2008; Mandelbaum et alj|2 008). 

In this letter we study the total mass-density profile of 
massive early-type galaxies inside one effective radius, us- 
ing the full SLACS sample of 58 gravitational lens systems 
with high-fidelity HST-ACS observations. We examine the 
intrinsic scatter in density slopes that is allowed by the sam- 
ple, whether the slope correlates with other global parame- 
ters and we set limits on the level of orbital anisotropy in 
these system. In Section 2, we present the logarithmic den- 
sity slopes of the SLACS early-type lens galaxies, based on 

^ www.slacs.org 



2 



Koopmans et al. 



two different methods, one based on lensing and dynamics 
and one based on scaHng relations, that explicitly include 
their average density profile as a free parameter. Comparing 
these two values allows us to set limits on their average orbital 
anisotropy. In section 3, we summarize our results and con- 
clusions. Throughout this paper, we make use of the samp le of 
58 SLACS single-lens systems from lBolton et al.l ( l2008ah and 
take all quantities from that paper. If not mentioned otherwise, 
all masses are in units of 1O"^M0 . We assume Q.^ - 0.3, 
Qa = 0.7 and Hq = 100 /i km/s/Mpc. 

2. THE DENSITY PROFILE OF MASSIVE ETGS 

To determine the logarithmic slope of the total density pro- 
file, we use two alternative methods: (i) through combining 
SDSS-based stellar velocity dispersions and lensing-based to- 
tal masses, and (ii) through scaling relations between lumi- 
nosity, mass and rescaled radius, which does not require a 
measured stellar velocity dispersion. Whereas a single-power 
law m as s model is an a pproximation, Koopmans & Tre^ 
(I2003h : iTreu & Koopmansl (120041) : iBarnabe et alJ (120091) 
show that two component models can robustly be approxi- 
mated by a single power-law component given current data 
quality. 

2. 1 . Derivation from Lensing and Stellar Dynamics 

T o derive th e logarithmic density slopes, w e fol - 
low iTreu & Ko opmans (2002); Koopmans & Tre^ (l2003h : 
iTreu & Koopmans (2004); Koopm ans et al.l (2006). First the 
lensing mass in side the critical curve s is determined from the 
lens models in iBolton et al.l ( l2008ah . whic h is nearly invari- 
ant under changes in the density profile (iKochanekl 1 1 99 Ih . 
hence the assumed density profile during th i s step is not rel- 
evant in further steps (see 'Koop mans et al.l (l2006l) ). Subse- 
quently we solve the spherical Jeans equations - assuming 
this mass as external constraint and an Einstein radius equal 
to the deflection angle of an equivalent spherical mass distri- 
bution - for a lum in osity-densit y profile that follows either a 
iHernquistI (Il990h or lJafi^d (IT983 ) profile, embedded as trace- 
component inside the total density profile ptot <^ ''^c (i.e. 
a power-law profile). The half-light radius of the projected 
luminosity -density pr ofile is set equal to the observed effec- 
tive radius (iBolton et al. 2008a). We take seeing (FWHM=1.5 
arc sec) into account. We vary the slope y'^^^ over a range 
of 1.1 to 2.9 and compare the predicted velocity dispersion 
inside the 3-arcsec diameter SDSS fiber with the observed 
value. The error in the measured velocity dispersion is by 
far the most dominant source of uncertainty. Hence, the like- 
lihood dPjdy'^^ oc e^^'l^ is determined from the;^'^ mismatch 
between the model and observed velocity-dispersion values. 

The results in Figure [T] show that most values lie around 
a slope of two, which is that of an isothermal mass profile 
(p oc r"^). A joint analysis of the sample yields: 

(tLd) = -«'log(ptot)/«'log(r) = 2.085!°:°25 (gg^^ ^L), (1) 

in the range of radii 0.2 - 1.3 /?eff and for /?,. = 0. The depen- 
dence on orbital anisotropy is small with the slopes varying 
mildly over y6,- = +0.50 (see Fig.[T]i. Based on changing the 
luminosity density profile, seeing, etc., we estimate an ~0.1 
systematic error. We note two points: (i) a more detailed two- 
dimensional kinematic analyses of six SLACS galaxies (see 
ICzoske et al.ll2008l: IBarnabe et all l2009) agree with these re- 
sults, and (ii) co mparing the densi t y slop es of the 14 systems 
that overlap with iKoopmans et al] (l2006l) . we find an average 
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Fig. 1. — The logarithmic density slopes of 58 SLACS early-type galaxies 
(thin solid curves). The filled red curve is the joint likelihood of the ensemble- 
average density slope. The histogram indicates the distribution of median 
values and the dotted Gaussian curve indicates the intrinsic scatter in y[ ^ 
(see text for details). We assume a Hemquist luminosity-density profile. Tne 
small dashes indicate the shift in the ensemble-average density slope for jS^ = 
+0.50, +0.25, -0.50 - 0.25 (left to right), respectively. Note the reversal of 
the fir = -0.50 and -0.25 dashes. The vertical solid line and gray region 
indicates the best-fit value and 68% CL interval, respectively, of the average 
density derived from scaling relations. 

increase of y'^j^ by 6%. This difference can be attributed to 
minor model improvements, the use of better HST images, 
leading to an average decrease of /?etf by 13%, and an im- 
proved derivation of the stellar velocity dispersion, leading to 
an increase by +3%. In particular the latter leads to an aver- 
age increase in y[^jj, explaining most of the difference. Hence 
currently we are limited by systematics. 
An intrinsic spread of o-y - 0.20^qq2 (i.e. 

o"y'/<rLD> = O-lO-SSop 68% CL) is derived, assuming Gaus- 
sian i ntrinsic and error distributions (see iKoopmans et al.l 
'2006'; 'Barnabe et al. "2009) consistent with the scatter found 
in Koopmans et al. (2006); Jiang & Kochanek (2007); 
Barna be et al.l (120091) . Despite the uniformity of the sample, 
differences between ga laxies are present, wh ich could partly 
be physical (see e.g. iGerhard et alj 1200 lb , partly due to 
systematics, or due to small uncorrected contributions from 
the environmen t and large scale structure (Auger 20081 
ITreu et alJl2009l: iGuimaraes & Sodre 2009). Conservatively 
it should therefore be regarded as an upper limit on physical 
variations. 

2.2. Derivation from Scaling Relations 

A second method to derive the ensemble-average density 
profile is to assume a scaling relation between the observ- 
ables, luminosity, effective radius and Einstein radius and Ein- 
stein mass: 

Meinst ( 



a log(Leff) = log 



eif 



+ 6 



(2) 



This relation assumes that the density profile scales as a 
power-law with p oc r"''sR and that mass and light scale with 
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Meff oc L«^, with Meff = Meinst (/?eft7/?einst)^^ The idea is 
that for a fixed luminosity, M^nsi scales with ^einst/^eff and 

thus provides an ensemble-averaged density slope^ 

Thi s approach is slightly different from iBolton et alj 
(|2008a), who assumed a SIE or const ant M/L mass profile 
in deriving the MFP, but similar to R usin & Kochanel^ (l2005l) 
who focused on the FP and density slope. It allows us to de- 
rive the density slope independently fmm assumptions about 
stellar dynamics. The resulting values are a = 1.363 + 0.056, 
6 = -0.959 + 0.050 and for the best-fit ensemble-average log- 
arithmic density slope 

(rsR>= 1.959 + 0.077. (3) 

This result is close to that derived based on dynamics models, 
although it assumes nothing about the dynamical structure of 
these galaxies (e.g. isot ropy). The value of a = l/rj', with rj' 
as defined and given in lBolton et al.l (12008 a). The diff'erence 
between {~y'gJ^) and ("Xl]-,) of ~ 0.1 irnplies that on average 
anisotropy can not be very large (see ^2Ai . 

2.3. Correlations of Slope and Galaxy Properties 

To assess whether y'^^^ correlates with global galaxy quan- 
tities or cosmic time, we plot them against effective mass, 
rescaled (/?einst/^eff) radius and redshift. The results are 
shown in Figure 2. We find the following linear gradients: 

dy'/dcTsm^(-l-6 + 83)-lO-^ 



dylJdM^s = (-0.6 + 14.5) • lO^^ 
dy[ Jd(R,inJRes) = (0.10 + 0.16) 



(4) 



[ dylJdz = (-0.05 + 0.43), 

with cTsiE in units of km/s and Mejf in units of 10'° M©. No 
weights on the points are used to avoid the brighter low- 
redshift galaxies from dominating the fits. We find no cor- 
rel ations at any sign ificant level. Similar res ults were foun d 
by iKooprnan s et al.' (200 6J and recen tly by 'A uged (l2008l) : 
iTreu et all (12009): Guima ries & Sodrel 12009). based on sfiid- 
ies of their environment and mass along the line-of-sight. 
Overall, the inner regions {R < /?eif) of massive early- 
type (> Lt) galaxies are remarkably homologous and sim- 
ple. This result is unlike some of the more detailed results 
found by the SAURON collaboration for < galaxies (e.g. 
ICappellarietan 120061: lEmsellem et"ani20()7h . which show a 
range of different kinematic signatures (e.g. counter-rotating 
cores, triaxiality, fast versus slow rotators, etc.). We find that 
these complexities, at least at the high-mass end, does not 
seem to affect the derivati on of their mass dis tributions and 
scaling relations (see e.g. iBarnabe et al.l 120091 for a discus- 
sion), although the overlap between the sampl es in mass is 
small. The analyses of six SLACS lens system (Czosk e et al.l 
12008: Barnabe et al. 2009). based on VLT VIMOS-IFU data, 
shows that the majority of the SLACS galaxies are slow rota- 
tors predominantly supported by random stellar motions. 

2.4. Limits on Orbital Anisotropy 

Of the two methods to derive the average density slope, one 
depends on dynamics and the other does not. In the former, 
we assumed fir - 0. If the galaxies, however, on average have 
a different value of y6,-, the average value of y'^^ can be either 
too large or to small, since we would attribute a higher/lower 
dispersion to a steeper/shallower density slope and not to ra- 
dial/tangential orbital anisotropy. 
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Fig. 2. — The median values of y' versu.s global galaxy quantities and 
redshift. The blue symbols are SO galaxies and the magenta colored system 
is an E/SO galaxy. The dashed curves is given at y'^j^ = 2.085 for reference 
with the gray box being ±10% intrinsic scatter. The thin red line is the best 
linear fit (i.e. shown curved in the log-plots). 

The independence of the second method from dynamics al- 
lows us therefore to place constraints on the average orbital 
anisotropy of the ensemble of galaxies (/J^). As indicated in 
Fig.l, good agreement is found for a positive value of (/?r) ~ 
0.45 ± 0.25, in agreement with results from iGerhard et al.l 
(1200 ll) . The galaxies are at most r nildly radial anisotrop ic, as 
also found in lCzoske etal] (l2008l) : lBarnabe et al.1 (120091) from 
a more detailed two-integral and two-dimensional analysis of 
SLACS systems. 

3. RESULTS & CONCLUSIONS 

We have presented a complete analysis of the inner mass 
density profile of 58 early-type galaxies from the SLACS sur- 
vey, focusing on their logarithmic density slope. We find the 
follow results for galaxies with Meif > 3 • IO'^Mq: 

• Based on lensing and stellar dynamics constraints, the 
average inner logarithmic density slope in their inner 
0.2-1.3 Reft is (Tld) = 2.085!°;°25 (random error on 
mean; 68% CL; for isotropic orbits, Pr — 0), +0.1 
(syst.) and a small (<0.1) dependence on anisotropy for 
a range of /?,• = ±0.50. An intrinsic scatter in y'^^, be- 
tween individual galaxies, of cry < 0.20!j;{;^ (68% CL) 
is still allowed by our data. This should be regarded as 
an upper limit on physical variations, since it must in- 
clude residual systematic effects. However, we believe 
part of it to be due to real di fferences between E TG 
density profile slopes (see also lBarnabe et al.ll2009l) . 

• A second independently-derived value of the average 
inner logarithmic density slope is found from scaling 
relations, yielding {y'^^d ~ 1-959 ± Q.Oni . This value is 
completely independent from dynamics and therefore a 
robust sanity check of the lensing plus dynamics results. 

• Since {y'^^ is independent of orbital anisotropy, we can 
set a weak limit (/?r) = 0.45 ± 0.25 on the average 



4 



Koopmans et al. 



anisotropy of these systems. This shows that massive 
galaxies are at most mildly radial anisotropic. 

• No correlation of y[^p is found with either galaxy mass 
or redshift, nor with radius over which this slope is mea- 
sured (0.2-1.3 Rett), implying that these results are gen- 
uine and widely applicable. 

• Taking the density slope into account, an increase in lu- 
minosity is found with increa sing galaxy mass: M etr 
^i363±{).056 agreement with lBolton et all ( l2008bl) . 

Based on these numerical results, we conclude that mas- 
sive early-type galaxies with total masses Meff > 3 ■ 10"^Mq 
are structurally and dynamically very similar in their inner re- 
gions (roughly one effective radius), over a look-back time 
of about 4 Gyr, with stellar and dark-matter adding up to a 
combined close-to isothermal density profile but having some, 
although relatively little, intrinsic scatter between their log- 
arithmic density slopes. This bulge-halo conspiracy occurs 
despite (i) a clear increase in their total M/L inside one ef- 
fe ctive radius w ith galaxy mass (as found in this paper and 
in iBolton et al.i r2008b). leading to a tilt in the FP, and (ii) a 
very complex hierarchical formation history. Moreover, from 
the agreement between two different determinations of their 



average density slope - one dependent on dynamics and the 
other not - we also find that these galaxies are also at most 
mildly radially anisotropy. Our results are based on single- 
component mass models, which provide a good description 
of the available data for systems. In forthcoming publica- 
tions, we will address two-component models where stellar 
and dark matter are modeled separately. 

Our results on the near homology, isothermality and 
isotropy of massive early-type galaxies provide new con- 
straints on theoretical models and numerical simulations of 
the formation of early-type galaxies, their subsequent evolu- 
tion and the understanding of the FP. Whereas these mod- 
els/simulations should match these scaling relations, they 
should also match their intrinsic scatter 
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